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The magneto-thermal conductivity k is calculated for the vortex state of UPd2Al3 by assuming 
horizontal gap nodes. The Green's function method we employed takes into account the effects of 
supercurrent flow and Andreev scattering on the quasiparticles due to Abrikosov's vortex lattice 
order parameter. The calculated angular dependence of k vv for field rotation Go in the ac-plane 
depends strongly on field strength H , impurity scattering, anisotropy of the Fermi velocity, and 
temperature. For finite temperatures and the clean unitary scattering limit we get qualitative 
agreement with recent experiments for all four proposed gap functions having horizontal line nodes 
at ck z = , ±7r/4 , and ±7r/2 . 



Angle-dependent magneto-thermal conductivity is a powerful tool for determining the nodal structure of the gap 
in unconventional superconductors. Recently, the thermal conductivity has been measured in the heavy-Fermion 
superconductor UPd2Al3 for a variety of magnetic field orientations. 1 The thermal conductivity K yy displays two- fold 
oscillations when the magnetic field H is rotated in the ac-plane, while no oscillations are observed when H is rotated 
within the basal ab-plane. These results provide strong evidence that the gap function A(fe) has horizontal line nodes 
orthogonal to the c-axis. Four gap models have been proposed and denoted in Ref. 1 as types I - IV, respectively: 
A(fe) oc sinx, cosx, sin 2%, and cos 2% (with \ = ck z )- Cooper pairing mediated by magnetic excitons yields the 
highest T c for the model cosx- 2 The magneto-thermal conductivity for these models has been calculated by the 
method of the Doppler shift of the energy of quasiparticles due to the circulating supercurrent flow of the vortices. 3,4 
This effect depends sensitively on the angle between H and the direction of the nodes of the gap. Comparison of 
these Doppler shift results with the data indicates that the model type IV, A(fe) = A cos 2\ , gives the most consistent 
description of the experiments of Ref. I. 3 A new Doppler shift approximation however yields the opposite result, i.e., 
that the models of types I - III are in better agreement with experiment. 4 

The purpose of the present paper is to calculate the magneto-thermal conductivity for horizontal gap nodes by 
another method which takes into account, beside the effect of the supercurrent flow, the Andreev scattering off the 
vortex cores. 5 These effects are calculated from the real and imaginary parts of the Andreev scattering self energy for 
the quasiparticle Green's function in the presence of Abrikosov's vortex lattice order parameter. This method was 
applied to Sr2Ru04 by assuming vertical and horizontal nodes of the superconducting gap. 6 An equivalent method 
based on the quasiclassical equations and linear response theory 7 provides much more compact expressions for the 
density of states and thermal conductivity. 8 We have shown that the latter expressions yield very nearly the same 
results as the original expressions derived in Rcfs. 5 and 6 (see Ref. 9). 

The Pesch-approximation 7 yields for the spatial average of the normalized density of states: 

N(lj, k) = Reg(u, k) ; g(u, k) = {l + 8|A(fc)| 2 [A/v ± (k)] 2 [1 + z w(z)\y 1/2 (1) 

where 

z = 2[w + iS i (o;)][A/«x(fc)] ; A = (2eH)~ 1 ^ 2 . (2) 

Here v±(k) is the component of the Fermi velocity perpendicular to the magnetic field H , and Sj(w) is the self energy 
for impurity scattering which is calculated self-consistently in the t-matrix approximation. 9 The integrand of the ui- 
integral for the thermal conductivity is proportional to (y) cosh -2 (wp) Keg(u),k) times the transport scattering 
time r(w,fe) where 8 

-* =RaEi ( t ,) + 2^[AMW]lA W |» ggM^^ ■ (3) 
2r(w,0) Reg{uj,k) 

The first term in Eq.[3] is the scattering rate due to impurities which reduces in the normal state to F = l/2r n . The 
second term is the Andreev scattering rate due to the vortex cores. 
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For the heavy-Fermion superconductor UPCI2AI3 we consider only the cylindrical Fermi surface and approximate 
the Fermi velocity by 



v F = v a [cos <f> 



e a + sin0e 6 + e sin2xe c ] , (x = ck z , -tt/2 < \ < +n/2) 



(4) 



where e = v c /v a . Furthermore we assume the following gap function with horizontal nodes: 



A(fc) = Acosx; (X = cfc z ); A = A [l - h 2 } 1/2 ; (h = H/H c2 ). 



(5) 



For a field direction given by polar and azimutal angles 9q and <j>o , i.e., 



H = H [sin 6*o cos <fio e a + sin O sin <j> + cos 9 e c ] , 



(6) 



the component of vp perpendicular to H becomes: 




sin O cos(0 — <fio) + £ cos #0 sin 2%] 
[1 + e 2 sin 2 2 X ] 



21 1/2 



V_L(fc) = V F 



(7) 



From this expression it is obvious that the combination A(k)/v±(k) occuring in Eqs. (1) and (3) depends sensitively 
on the direction of H relative to k and thus relative to the gap nodes. For this model we have calculated K yy /n n 
for field orientations H in the xz-planc (0o — 0; < 8q < ir). For the anisotropy parameter e = v c /v a in 
Eq.(4) we take the value e = 0.6 . We first consider an impurity scattering rate S = r/A = 0.024 and a scattering 
phase shift p = 0.9(7r/2) . In Fig. la we show n yy /n n versus the reduced field h for 9 a = (solid curve) and 
8q = ir/2 (dashed curve) for T = (u> = 0). We have also plotted the relative amplitude of the oscillation, 
Ak = [K yy (9o = 0) — Kyy(0o = tt/2)]/ K yy (9o = 0) . In Fig. lb we show two examples of the oscillation as functions 
of the polar angle 9q for the fixed field strengths h = 0.05 and h = 0.3 . Our results for the unitary limit (phase 
shift p = tt/2 ) are nearly the same as those shown in Fig. 1. Different results are obtained in Ref. 4 for phase shift 
p = 0.9(7r/2) where a small region (denoted by I in Rcf. 1) exists in which n yy decreases with increasing h and the 
amplitude of the oscillation An becomes negative. This would be in agreement with the data of Ref. 1. 

The measured angular variation of K yy (9o) is decomposed into a constant and a term with twofold symmetry with 
respect to the 9 rotation, i.e., C yy cos2# - 1 In the region denoted by II in Rcf. 1, where n yy depends linearly on H , 
C yy is positive. In the small region I near H — , where n yy decreases with increasing H , the measured amplitude 
Cyy becomes negative. 

We have not been able to confirm the results of Rcf. 4 where the origin of the sign change of An at low fields is 
the tiny deviation of the impurity scattering phase shift 0.9(7r/2) from the unitary limit 7r/2 . We show now that this 
sign change is actually due to the finite temperature. We estimate the finite temperature effect by considering finite 
values of O = uj/A where u occurs in Eqs. (2) and (3) and is the integration variable for the ^-integral yielding k . 
Since this integrand is strongly peaked at uo/T ~ 2.4, the temperature corresponding to f2 is approximately given 
by £1 ~ 2.4(T/A ) ~ T/T c . We have verified that this is a good approximation by numerical computation of the 
w-integral at finite T . As an example, in Fig. 2a we show n yy /K n versus h for f2 = 0.3 for the impurity scattering 
rate 5 = 0.024 and phase shift tt/2 . The solid curve refers to the angle 9 — and the dashed curve to 9 = tt/2 . 
One observes now a small range from h = to about h = 0.1 where the amplitude An of the oscillation is negative. 
In the enlarged scale of Fig. 2b it is seen more clearly that in this region K yy j n n bends upwards for decreasing h , 
which also occurs in region I in the experiments. In Fig. 2c we have plotted the oscillation of n yy as a function of the 
polar angle 9 for fixed field strengths h — 0.05 , 0.3, and =0.5. One sees that oscillations of approximately twofold 
symmetry in 9q occur. In the region of approximately linear field dependence of K yy (see Fig. 2a), which corresponds 
to region II in the experiments, the amplitude Ak corresponding to C yy in Ref. 1 is positive. In the region of low 
fields, corresponding to region I of Ref. 1, An is in fact negative for this finite temperature. It is important to note 
that the experiments are carried out at a finite temperature of about 0.4 K. 1 

We discuss now the dependence of our results on the various parameters. First we consider the effect of impurity 
scattering. For a larger scattering rate, for example, 5 = T/Ao = 0.1 , and the unitary limit, we obtain for Q = 0.3 
the functions K yy (h) and Ak shown in Fig. 3a. Comparison with the results shown in Fig. 2a for 5 = 0.024 , and all 
other parameters the same, shows that the slope of K yy in the linear range of region II is larger, while it is smaller in 
the region, denoted by III in Ref. 1, where K yy /K n rises steeply to one. This field dependence in Fig. 3a is in better 
agreement with experiment. However, no region I occurs where k vv decreases with h and Ak becomes negative. The 
function K yy (9o) exhibits an approximate twofold symmetry in 9 n for all values of h. For S = 0.1 we have to go to 
higher values of Q to obtain a sign change in An near h — . This is shown in Fig. 3b where Ak and k vv are plotted 
vs h in an enlarged scale for Q = 0.3, 0.4, and 0.5. One sees that for il = 0.4 (0.5) , Ak becomes negative below 
h ~ 0.06 (0.18) . These regions of negative Ak correspond to region I because K yy bends upwards for decreasing h . 
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Another important parameter is the phase shift for impurity scattering. We have already pointed out that, contrary 
to the results of Ref. 4, we obtain practically no difference between the results for the unitary phase shift limit ir/2 
and 0.9(7r/2) . However, for the Born limit, phase shift zero, the curve for K yy (h) is much flatter in the lower /i-range 
than the curves for the unitary limit in Figs. 2a and 3. Also, the amplitude Ak is negative and |Ak| is much larger 
over the entire range of h . These results in the Born limit are in agreement with those obtained in Ref. 4. 

We come now to the most crucial parameter, e = v c /v a in Eq.(4), which determines the anisotropy of the Fermi 
velocity for the corrugated cylindrical Fermi surface sheet. For the range of ^-values between about e = 0.5 and 0.7 
we obtain nearly the same results as those shown for e = 0.6 . These values agree approximately with the e value 
which is obtained according to Ref. 4 if we take for v 2 and v 2 the average values over the corrugated cylindrical 
Fermi surface where (v c /v a ) 2 = (l/2)(v c0 /v a0 ) 2 and a = (v c0 /v a0 ) 2 ~ 0.69 is the appropriate value for UPd 2 Al 3 . 
In addition to the order parameter A oc cos\ (type II), we also obtain, for e between 0.5 and 0.7, nearly the same 
results for the other order parameters proposed in Ref. 1, namely, A oc sin% (type I), sin 2% (type III), and cos 2% 
(type IV). In all cases \ — °kz ■ However, for Q small, as e increases from 0.7 to about 0.84, the minimum of K yy (0o) 
at the field angle 6q = (tt/2) changes to a smaller maximum and two minima develop at the intermediate angles of 
about 8o ~ ±1.3(7r/4) . It should be pointed out that, for the type IV gap, the minima already start to form just 
above e — 0.6 and become more pronounced than for the type I - III gaps. For increasing f2 , the small peak at n/2 
is reduced and becomes a flat, broad minimum at about = 0.3 . Furthermore, for f2 above 0.26, the amplitude Ak 
of the variation of n yy (9o) becomes negative in a small range near h = . This behavior is similar to that of Ak 
shown in Fig. 2 for e = 0.6 and Q = 0.3 and again we see the importance of finite Q (temperature) in reproducing 
the experimental results. We remark that, for SI = 0.3 and e = 0.6 , Ak is negative below h = 0.05 for the type IV 
gap while it is negative below h = 0.1 for the type II gap. 

We also find that, for ft larger than 0.4, a sign change of Ak occurs in a small range of h just below h = 1 . In 
this connection we remark that experiments yield a negative sign of the amplitude C yy in a region denoted by III 
in Ref. 1 where K yy rises steeply with H . This negative sign has however been attributed to the anisotropies of the 
Fermi velocity and the upper critical field H C 2 . Since we have neglected here the anisotropy of H C 2 it may be that 
our calculated values of |Ak| in region III are much smaller than the measured values. 

Let us compare our results with those obtained by the Dopplcr shift method. For the clean unitary scattering 
limit and intermediate field strengths we obtain for the angular variation of n yy {0o) a positive amplitude Ak = 
[K yy (9o = 0) — K yy (9 = n/2)]/K yy (6o = 0) for all four of the gap models of types I - IV proposed in Ref. 1. The first 
Dopplcr shift calculation 3 yielded a positive Ak for the gaps cos 2% and sin% and a negative Ak for cosx . Another, 
more recent, Doppler calculation 4 finds a positive Ak for sin% , cosx , and sin 2% whereas the variation for the gap 
cos 2% is completely different from the form cos2# and exhibits pronounced minima at the angles 9 = ±51° . It is 
interesting that we obtain similar structure in K yy (0o) with minima at intermediate angles 9 n ~ ±1.3(7r/4) for all four 
gap models I - IV if the anisotropy parameter e — v c /v a is increased from 0.6 to 0.84. It is not surprising that our 
analytical expressions yield different results than the Doppler shift method because we take the Andreev scattering 
of the quasiparticles by the vortex cores into accoount. It has been shown by comparison with the exact results 
of the quasiclassical Eilenberger equations (see Fig. 7 in Ref. 10) that the analytical expressions 5,7 provide good 
approximations over the whole field range from H C 2 to H c i wheras the Doppler shift method is a good approximation 
only at low fields. However we want to caution that in the low field regime the higher Fourier components of the 
Green's function with respect to the vortex lattice vectors 5 might lead to important corrections. Therefore a fully 
self-consistent numerical calculation, in particular for the internal field, would be necessary 11 in order to obtain more 
reliable results. This means that the method of Ref. 11 for the density of states should be extended to the calculation 
of the Andreev scattering from the supcrfluid flow in the single vortex regime which is important for the thermal 
conductivity k . The large effect of this Andreev scattering can be seen from the fact that, at T — and impurity 
scattering phase shift = 0.9 (tt/2), k increases for small fields while the Doppler shift method together with a scattering 
lifetime due to impurity scattering only yields an initial drop of k for small fields. 4 We stress again that only at finite 
temperature (ft) do we obtain an initial drop of k together with a sign reversal of the oscillation amplitude Ak as 
shown in Figs. 2b and 3b. 

In summary, in the clean unitary scattering limit, and for not too large anisotropy e of the Fermi velocity, we obtain 
angular variations of K yy (9o) of approximately twofold symmetry in the field rotation angle 9q for all of the four 
proposed gap models having horizontal line nodes at ck z = \ = , ±7r/4 , and ±ir/2 . The amplitude Ak is positive in 
the range of field values H where n yy is nearly linear va H . At finite temperature Ak becomes negative at small fields 
which is in agreement with experiment. The range of field values with negative Ak is smaller for the gap of type IV 
(cos 2%) than for types I - III. Our results depend crucially on the anisotropy parameter e = v c /v a for the corrugated 
cylindrical Fermi surface. For values of e above about 0.7 we find a small maximum at 9o = 7r/2 and two minima at 
about 6*o — ±1.3(7r/4) in the curve of K yy (0o) for the gap models I -III. For the gap model IV this critical value of e is 
smaller, about e — 0.6 . It should be pointed out that even for a constant gap the density of states depends strongly 
on the shape of the Fermi surface and the field orientation. 12 Thus we can say that the angular variations of n yy {9 ) 
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are determined by two competing effects: the positions of the gap nodes and the amount of anisotropy of the Fermi 
velocity given by the shape of the Fermi surface. Here the latter effect appears to predominate. 

In conclusion, our calculations indicate that at present one cannot determine which one of the proposed gap models 
with horizontal line nodes is most consistent with the measured magneto-thermal conductivity of UPd2Al 3 . The gap 
model cos \ however seems to be somewhat more favorable than cos 2% . Our calculations also indicate the importance 
of accounting for finite temperature and lead us to expect that the value of the field at which the amplitude of the 
oscillation An changes sign should decrease with decreasing temperature. It is thus desirable to have measurements 
available at lower temperatures in very clean samples. It will also be interesting to see whether additional structure, 
minima, for example, occur in the form of K yy (9o) ■ 
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Fig. 1a 
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FIG. 1. la) Thermal conductivity K yy for the gap model I vs reduced field h = H/H C 2 for field direction 
8q = (solid curve) and 60 = it/ 2 (dashed curve). Also shown is the relative amplitude of the oscillation 
Ak = [K yy (9o = 0) — K yy (9() — n/2)]/n yy (9o = 0) . The impurity scattering rate is 8 — F/Ao = 0.024, the phase shift is 
0.9(7r/2) , and the anisotropy parameter e = v c /v a = 0.6 . 
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Fig. 1b 
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FIG. 2. lb) Angular variation of k vv vs polar angle 6o for field strengths h = 0.05 and h = 0.3 . 
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Fig. 2c 
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FIG. 5. 2c) Angular variation of n yy (do) for field strengths h — 0.05 , 0.3, and 0.5. 
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Fig. 3a 




Fig. 3b 




